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Within a pQCD inspired kinetic parton cascade we simulate the space time evolution of 
gluons which are produced initially in a heavy ion collision at RHIC energy. The inelastic 
gluonic interactions gg <-> ggg do play an important role: For various initial conditions it is 
00 | found that thermalization and the close to ideal fluid dynamical behaviour sets in at very 
early times. Special emphasis is put on color glass condensate initial conditions and the 
'bottom up thermalization' scenario. Off-equilibrium 3^2 processes make up the very 
beginning of the evolution leading to an initial decrease in gluon number and a temporary 
avalanche of the gluon momentum distribution to higher transversal momenta. 

Ph! 1. INTRODUCTION 

Qh' 

It had been demonstrated that the measured momentum anisotropy parameter v 2 at 
RHIC energy can be well understood if the expanding quark-gluon matter is described 
by ideal hydrodynamics. This important finding suggests that a strongly interacting 
^ . and locally thermalized state of matter has been created which behaves almost like a 
perfect fluid. On the other hand, the initial situation of the quark-gluon system is far 
from thermal equilibrium. It is thus important to understand how and which microscopic 
partonic interactions can thermalize the system within a short timescale and can be 
responsible as well for its (nearly) ideal hydrodynamical behaviour. 

A traditional way to study thermalization of particles is to carry out microscopic trans- 
port simulations. A standard parton cascade analysis incorporating only elastic (and 
forward directed) 2^2 collisions described via one-gluon exchange, shows that thermal- 
ization and early (quasi-)hydrodynamical behaviour (for achieving sufficient elliptic flow) 
can not be built up or maintained, but only if a much higher, constant and isotropic cross 
section a e ff ~ 45 mb is being employed [H]. The gluons being treated as semi-classical 
degrees of freedom would then aquire a collional width T m na e ffV re i being much larger 
than the typical energy E « 3T. The gluons would then resemble very broad excitations 
and can not be considered at all as quasi-particles, if such large and constant cross sections 
are employed. 

In contrast, a kinetic parton cascade algorithm [[3] has beeb deneloped with pertur- 
bative QCD inspired processes including for the first time inelastic ('Bremsstrahlung') 
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Figure 1. Time evolution of nomalized 
transverse energy at midrapidity. Results 
are obtained from simulations with differ- 
ent initial conditions of gluons for a cen- 
tral Au Au collision. 
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Figure 2. Time evolution of the elliptic 
flow v 2 at midrapidity for various impact 
parameter b. 



collisions gg <-» ggg. The multiparticle back reation channel is treated fully consistently 
by respecting detailed balance within the same algorithm. The three-body gluonic inter- 
actions are described by the matrix element 

lM "" 1 = (— (qi + m^ J Ui[(^-q.) 2 + ^J ( " * ^ ( ) 

where an effective Landau- Pomeranchuk-Migdal suppression is considered. A g denotes 
the gluon mean free path, which is given by the inverse of the total gluon collision rate 
A g = 1/Rg. This leads to a lower cutoff of k± and to an effective increase of the col- 
lision angles. Typical Debye screened total pQCD cross section scale roughly like the 
inverse temperature squared, ~ ^/T 2 , so that the characteristic ratio Y/E should be 
(significantly) smaller than 1, as long as the coupling stays small enough. The possible 
importance of incorporating inelastic reactions on overall thermalization was raised in 
the so called 'bottom up thermalization' picture [Ej- It is intuitively clear that gluon 
multiplication should not only lead to chemical equilibration, but also should lead to a 
faster kinetic equilibration. This has been demonstrated in detail in 

A first picture of the production of the primary partons at the very onset of a heavy 
ion collision is based on a free superposition of minijets being liberated in the individual 
semihard nucleon-nucleon interactions with transverse momentum being greater than a 
certain cutoff Po [El- Conditions of a color glass condensate fot the initial partons will 
also be given and further detailed in the next section. 

Phenomenologically the cutoff po can be chosen in a way to fit the (final) dEr/dy 
as seen in experiment. A value of po = 1.4 GeV would meets a final dE^/dy ~ 625 
GeV. Thus the minijets production at po = 1.4 GeV seems to give an appropriate initial 
condition of gluons at RHIC. In Fig. ^the transverse energy at midrapidity is depicted, 
obtained from simulations with the different initial conditions, being normalized to the 
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Figure 3. The evolution of the local mo- Figure 4. Respective transverse momen- 
mentum occupation within a typical color turn spectrum in the central region, 
glass condensate initial condition taken at 
t = 0.4 fm/ c is shown for four subsequent 
times. 

energies at t = 0.5 fm/c. One sees a rather unique behavior of a decreasing, normalized 
dET/dy\ y= o, especially for initial conditions employing minijets with p = 1.3 — 1.5 GeV 
and CGC with Q s = 1.0 GeV. Comparisons to the normalized transverse energy per unit 
rapidity for an ideal hydrodynamical expansion show that the collective expansion due 
to the pQCD interactions is quasi ideal at 0.5 — 1.5 fm/c. At later times the decrease of 
dET/dy\ y= o slows down, since the collision rates become smaller, especially in the outer, 
transversally expanding region. 

Taking po = 1.4 GeV for the initial minijets, the parton evolution for noncentral col- 
lisions at RHIC energy is simulated in order to calculate the elliptic flow parameter t>2. 
Figure 121 shows the time evolution of v 2 extracted at midrapidity for various impact pa- 
rameter b. These calculations are still preliminary. The results gives strong indication 
that an early pressure is being built up within that pQCD inspired description. The sym- 
bols in Fig. |21 mark the time from which the energy density in the central region decreases 
below 1 Gev/fm 3 . If we take the t>2 values at that marked times as the contribution from 
the partonic phase, they lie well in the region covered by the experimental data. 



2. Thermalization of the color glass condensate 

For the initial gluon distribution from the saturation picture we employ an idealized 
and boost-invariant form [Ej and express it as 

a s l\ c Tf 

which is described by Q s , the momentum scale at which gluon distribution saturates. 
The boost-invariance leads to the equality of momentum and space-time rapidity, i.e., 
77 = y, for the initial gluons. The value of parameters are taken from [ 5j: N c = 3 for 
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Figure 5. Time evolution of particle num- Figure 6. Time evolution of the effective 
ber per rapidity for various saturation mo- temperature T = E/3N. 
menta Q s . 

SU(3), c = 1.3, a s = 0.3, and the corresponding formation time, at which the gluon 
distribution becomes dilute enough, is given as r = c/(a s N c Q s ). Hence, the cascade 
operates from r = 0.4fm/c for Q s = 1 GeV and r = 0.18fm/c for Q s = 2 GeV and 
Q s = 3 GeV. For the calculation depicted in Figs. El and |U the gluons are produced 
within a transverse radius of 6 fm (Au nucleus) and within \n\ < 3 longitudinally. Fig. 
E] states a first dynamical realization of the so called 'bottom up thermalization' scenario 
as advocated in [ The evolution of the momentum occupation is shown for four 
subsequent times sampled within a space-time rapidity interval of An = 0.1 and a central 
transverse region of R < 1.5 fm. One recognizes the population of the 'soft' gluons and 
a subsequent degradation of the 'hard' initial gluons. In addition, also harder gluons 
are produced. All this happens roughly within the first 1 fm/c. In the last picture we 
also see that all particles are clearly more centered around the origin, demonstrating the 
ongoing cooling and quasi hydrodynamical behaviour from 2 to 4 fm/c. In Fig. 0] the 
transverse momentum spectrum is given for various times. Comparing the spectra 
between 2.0 and 4.0 fm/c, the local system at the central region clearly is at equilibrium 
at that later times, and the spectrum steepens continuously, since the thermo dynamical 
work in outward direction cools down the system. 

Thermalization somehow resembles the idealistic 'bottom up scenario' [ |3] with low 
momentum gluons being produced and populated. On the other hand, a couple of im- 
portant differences do show up. In Fig. |5] the time evolution of particle number per 
rapidity for various saturation momenta Q s is shown. In these simulations strict Bjorken 
geometry in longitudinal direction is preserved by reflection at cylindrical boundaries [ 
H]. The number of gluons is initially slightly decreasing, although a strong parametric 
enhancement has been advocated in [ |3] due to Bremsstrahlung production, which is not 
observed in the present calculations. This is true for all chosen saturation scales. The 
initial gluons are seemingly oversaturated, if the system thermalizes rather immediately [ 
H] . Only if thermalization would happen on a much longer time scale, the gluon number 
would become undersaturated at a subsequent time, at which net gluon production would 
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Figure 7. Time evolution of the rescaled temperature T-tz and the momentum anisotropy. 



then take place. 

Also it is found [0] that indeed the ggg — > gg processes dominate the dynamics from 
the very beginning up to approximately 2 fm/c. This causes, in particular, a production 
of gluons with momenta much larger than the saturation scale Q s . This can also be seen 
from Figs. El and 01 where gluons with such high momentum do show up. This behaviour 
resembles an avalanche from low momentum to high momentum, which has recently also 
be found in classical Yang- Mills simulation [ 6j. Both phenomena might be connected and 
clearly deserve further investigation. It also seems that the gluon distribution at higher 
momenta above the saturation scale evolves more faster to an exponential shape than the 
low momentum gluons [Ej. The name bottom-up has then to be reversed. 

As a last point the extraction of a possible thermalization scale is briefly discussed [0]. 
Fig. IH1 gives the time evolution of the effective temperature T = E/3N for three different 
saturation scales. At the very beginning the temperature increases by a small amount 
because of the decrease in gluon number, but then decreases more or less in accordance 
with ideal 1-dimensional hydrodynamical expansion. This can be seen more explicitly 
in the upper Fig. [7| (Quasi-)Ideal hydrodynamical behaviour sets in at a timescale 
9 = 3.3 fm/c{Q s = 1 GeV), 9 = 1.9 fm/c{Q s = 2 GeV) and 9 = 1.2 fm/c{Q s = 3 GeV). 
Thermalization should set in slightly earlier. Kinetic equilibration is characterized by the 
time evolution of the momentum anisotropy, 2 < p 2 z > / < p\ > , which maintains a 
value of one at equilibrium. The momenta of gluons are initially directed transversely in 
local discs according to Eq. (J2J). Therefore the momentum anisotropy has a vanishing 
initial value. An exponential fit to the time evolution gives 9 = 2.5 fm/c(Qs = 1 GeV), 
9 = 1.5 fm/c(Q s = 2 GeV) (cf lower Fig. GJ, = 1.0 fm/c(Q s = 3 GeV), although 
equilibration happens faster than exponetial at the beginning. 
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